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Inverse Compton e pair cascade model for the gamma-ray 
production in massive binary LSI +61° 303 
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ABSTRACT 

We apply an inverse Compton pair cascade model for 7-ray production in the 
massive binary system LSI +61° 303 assuming that electrons are accelerated already 
inside the inner part of the jet launched by the compact object. 7-ray spectra, affected 
by the cascade process, and lower energy spectra, from the synchrotron cooling of the 
highest energy electrons in the jet, are calculated as a function of the phase of this 
binary system. 7-ray spectra expected in such model have different shape than those 
ones produced by electrons in the jet directly to observer. Moreover, the model predicts 
clear anti-correlation between 7-ray fluxes in the GeV (1-10 GeV) and TeV (> 200 
GeV) energy ranges with the peak of the TeV emission at the phase ~0.5 (the peak 
half width ranges between the phases ~0.4-0.9 for the inclination of the binary system 
equal to 60°, and ~0. 4-0.1 for 30°). The fine features of TeV 7-ray emission (fluxes 
and spectral shapes) as a function of the phase of the binary system are consistent 
with recent observations reported by the MAGIC collaboration. Future simultaneous 
observations in the GeV energies (by the GLAST and AGILE telescopes) and in the 
TeV energies (by the MAGIC and VERITAS telescopes) should test other predictions 
of the considered model supporting or disproving the hypothesis of acceleration of 
electrons already in the inner part of the microquasar jets. 
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1 INTRODUCTION 

The binary system LSI +61° 303 contains a massive Be type 
star, with a radius of r*=13.4 R0 and a surface tempera- 
ture T s = 2.8 x 10 4 K, and a compact object on a close orbit 
characterised by a semi-major axis a — 5.3?%, an ellipticity 
e = 0.72, and an azimuthal angle of the observer in respect 
to the periastron passage u> — 70° (Casares et al. 2005). The 
inclination angle of the binary system is not well constrained 
by the observations and ranges from 25° < i < 60° for the 
case of a neutron star and 8 < 25° for a black hole. In the 
radio wavelengths, the observed asymmetric extended struc- 
ture is usually interpreted as due to relativistic radio jets, 
with an estimated speed ~0.6c (Massi et al. 2004). This ra- 
dio source has been related by Gregory & Taylor (1978) to 
the COS B 7-ray source CG135+01 (Hermsen et al. 1977). 
CG135+01 has been also detected by the EGRET instru- 
ment above 100 MeV (3EG J0241+6103, Hartman et al. 
1999) with a hard spectrum (with a photon spectral index 
2.05 ± 0.06, Kniffen et al. 1997), and by the COMPTEL de- 
tector in the energy range ~0.75 to 30 MeV (with a photon 
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spectral index 1.95_ ' 3 , van Dijk et al. 1996). The analysis of 
different EGRET observations shows evidences of variabil- 
ity (Tavani et al. 1998, confirmed by Wallace et al. 2000) 
presenting a modulation with timescales similar to the or- 
bital period of LS I +61° 303 and the emission maxima likely 
near the periastron passage and phase 0.5 (Massi 2004). The 
Whipple group (Hall et al. 2003, Fegan et al. 2005) set an 
upper limit on the TeV flux from this source, 0.88 x 10" 11 
cm" 2 s" 1 above 500 GeV and 1.7 x 10" 11 cm" 2 s" 1 above 
350 GeV, which is clearly below an extrapolation of the 
EGRET spectrum. Very recently, the MAGIC telescope ob- 
served TeV 7-ray emission at energies above ~200 GeV with 
the power law spectrum (spectral index -2.6) and extend- 
ing up to a few TeV (Albert et al. 2006). This emission is 
observed after the periastron passage of the compact ob- 
ject. 7-ray emission from LSI +61° 303 has been interpreted 
in terms of the inverse Compton (IC) scattering model in 
which stellar radiation is scattered by electrons accelerated 
on the shock wave either created in collisions of the stellar 
and pulsar winds (e.g. Maraschi & Treves 1981), or by elec- 
trons accelerated in the jet and scattering of synchrotron 
and stellar radiation (e.g. Bosch- Ramon & Paredes 2004). 
The alternative hadronic model for this source has been 
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also discussed by Romero, Christiansen & Orellana (2005). 
The variable GeV gamma-ray emission reported from LSI 
+61° 303 by EGRET can be naturally explained provided 
that the production region of 7-rays lays inside the binary 
system, i.e. not far away from the massive star. Then, the 
optical depths for accelerated electrons (and produced 7- 
rays) are large enough for developing the IC e pair cas- 
cade in the radiation field of the massive star (Dubus 2006, 
Bednarek 2006 (B06)). For the geometry of the considered 
binary system, the radiation field is strongly anisotropic if 
considered with respect to the injection place of the elec- 
trons. Such general type of anisotropic IC e pair cascades 
in the radiation of the massive stars have been considered 
in the past by Bednarek (1997, 2000) and Sierpowska & 
Bednarek (2005), and recently applied to the binary systems 
LS 5039 and LSI +61° 303 by Bednarek (B06). Aharonian 
et al. (2006) also consider anisotropic cascades with the ap- 
plication to LS 5039. The basic features of such cascades 
differ significantly from the well-known IC e pair cascades 
in the isotropic radiation fields studied in detail in the past 
(e.g. Protheroe 1986, Honda 1989, Protheroe & Stanev 1993, 
Biller 1995). 

In our previous paper (B06), we considered 7-ray pro- 
duction in the cascade model assuming that electrons are in- 
jected isotropically at some distance from the massive star. 
Based on those calculations, we were able to predict at which 
orbital phases the largest fluxes of TeV 7-rays should be 
expected. In fact, these predictions have been successfully 
tested by the MAGIC observations (Albert et al. 2006). In 
this paper, based on the anisotropic cascade calculations 
(B06), we discuss a more specific model for the 7-ray pro- 
duction in LSI +61° 303. Here we assume that electrons 
are accelerated on shocks waves propagating along the jet 
launched from the vicinity of the compact object. They are 
injected along the jet with specific injection efficiency which 
does not depend on the phase of the binary system. We do 
not consider the effects of possible changes of the accelera- 
tion efficiency related to, e.g. variations in the accretion rate 
onto the compact object and neglect the possible radiation 
from the equatorial disk around the massive Be star in re- 
spect to its emission from the surface. In the present code we 
include the synchrotron losses of electrons in the jet which 
were not considered in previous work (B06). 

2 IC E ± PAIR CASCADE MODEL 

In the simple model considered in this paper it is assumed 
that electrons are accelerated continuously along the jet 
propagating from the compact object perpendicularly to the 
plane of the binary system. We can apply a simple scaling of 
the magnetic field along the jet assuming its conical struc- 
ture with the fixed opening angle a, 

B(x) = B dV /(l + az) « B dV /(az) for z > 1/q, (1) 

where a is assumed of the order of 0.1 rad, the scaling factor 
r\ < 1, and z is the distance along the jet from its base. We 
assume B d Ri3x 10 5 G, i.e. the order of magnitude applied 
also in other papers (e.g. Bosch- Ramon et al. 2006). The 
acceleration rate of electrons (with the Lorentz factors 7) in 
the jet can be parametrised by, 

(7) =^E/r L « 10 13 <eB eVs" 1 , (2) 



where E is the energy of accelerated electrons, £ is the ac- 
celeration efficiency, tl = E/eB is their Larmor radius, B is 
the magnetic field at the acceleration site (in Gauss) , e is the 
electron charge, and c is the velocity of light. The basic role 
in the above formula is played by the acceleration efficiency 
of particles £. The theory of shock acceleration process esti- 
mates the value of £ as ~ 0.1(v/c) 2 , i.e. ~ 0.04 for v = 0.6c 
(Malkov & Drury 2001). On the other hand, the presence of 
relativistic leptons with energies ~ 10 15 eV inside the Crab 
Nebula suggests that £ should be not very far from unity for 
relativistic winds. Therefore, we consider values of £ in the 
range 0.3 - 0.03. 

During the acceleration process, electrons lose energy 
mainly by the synchrotron and the inverse Compton (IC) 
processes. We are interested in the situations in which elec- 
trons can be accelerated to TeV energies. This is not easy to 
obtain relatively close to the base of the jet where the mag- 
netic field is large and the synchrotron losses very strong. 
Therefore, in the case of the IC process we consider dis- 
tances from the accretion disk for which the energy losses 
on the massive star radiation dominate and the IC losses 
on the disk radiation and the synchrotron photons can be 
neglected. 

Let us estimate the maximum energies of electrons in 
such shock acceleration process allowed by the energy losses. 
The energy loss rate by the synchrotron process is given by, 

A y „(7) = ^ttcttcUb"/ 2 « 2.07 x 10" 3 B 2 7 2 eV s~\ (3) 

where Ub = B 2 /8tt is the energy density of the magnetic 
field, 7 is the Lorentz factor of electrons, and or is the 
Thomson cross section. Electrons lose energy on the IC pro- 
cess in the Thomson (T) regime, if their Lorentz factors are, 

7 < 7T/kN = rnc 2 /3k B T ^2x 10 5 /T 4 , (4) 

(where T = 10 4 T4 K is the surface temperature of the mas- 
sive star which determines the black body spectrum, and fca 
is the Boltzman constant), and in the Klein-Nishina (KN) 
regime, for the Lorentz factors 7 S> 7t/kn- The energy loss 
rate in the Thomson regime is given by, 

Ac(t) = |iraTcI/ rBd 7 2 « 3.8T 4 V/r 2 eV s _1 , (5) 

where C/ rad = 4.7 x lQ li T%/r 2 eV cm -3 , and r is the distance 
to the centre of the massive star in units of its radius r*. The 
energy loss rate in the KN regime depends only logarithmi- 
cally on the Lorentz factor of electrons. We approximate 
these losses by (see e.g. Blumenthal & Gould 1970), 

P™{l) « Pic(7t/kn) In (4fc B T 7 /m C 2 - 2). (6) 

The maximum energies of accelerated electrons are deter- 
mined by the balance between the acceleration mechanism 
and by the most efficient mechanism for the energy losses. 
The comparison of Eqs. [5] and [3] gives the maximum allowed 
Lorentz factors of electrons due to the synchrotron energy 
losses, 

7s yn « 6 x 10 7 (£/B) 1/2 . (7) 

The maximum energies allowed by IC process in the Thom- 
son regime are (from Eq. |2Jand|SJl, 

7i T c w 1.6 x 10 6 (£B) 1/2 r/T 4 2 , (8) 
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Figure 1. The maximum Lorcntz factors of electrons as a func- 
tion of the distance, z, from the base of the jet (measured in units 
of the stellar radius), for different acceleration efficiencies: £ = 0.3 
(full curves), £ = 0.1 (dotted), and 0.03 (dashed). Specific figures 
show the results for two ratios of the magnetic field to accretion 
disk radiation energy densities counted at the base of the jet, 
equal to r\ = 0.1 (figures a and c), and 0.01 (b and d), for the 
periastron passage of the compact object (a and b) and for the 
apastron passage (c and d). These maximum Lorentz factors are 
determined by different energy loss processes. The comparison of 
the acceleration rate with the energy loss rate on the synchrotron 
process gives the straight parts of curves with dependence <x z 1 / 2 
(marked by syn). At some distances acceleration of electrons is 
saturated by energy losses on the inverse Compton scattering with 
stellar photons in the Thomson and the Klein-Nishina regimes 
(parts of curves corresponding to broad parabola dips, marked 
by IC). 

and in the KN regime are (from Eq. |5]and ^J, 

7l K c N n 1.5 x 10 5 (2 + e^ Br ' 2 ^)/T 4 . (9) 

Note that the limit given by Eq.|SJis only valid provided that 
Tic < 7t/kni i- e when the parameters of the jet and the 
massive star fulfil the following condition, T&/r > 8(^B) ' 5 
(from Eq. 0] and Ej. In any other case, the limit given by 
Eq.[!5]has to be taken. 

The maximum Lorentz factors, to which electrons can 
be accelerated at a specific part of the jet, are determined 
by the acceleration efficiency and the energy losses through 
synchrotron and IC processes. These energy loss processes 
depend on the conditions in the jet (magnetic field strength, 
acceleration efficiency), on the parameters of the massive 
star (its radius, surface temperature), and on the distance 
of the acceleration region from the massive star. In Fig. Q 
we plot the maximum Lorentz factors of electrons as a func- 
tion of the distance from the base of the jet (given by either 
Eq.Q or Eq.|S] or Eq.|5J, for selected values of the accelera- 
tion efficiency, £, and the ratio, r\, at the periastron and the 
apastron passages of the compact object on its orbit around 
the massive star. It is clear that the maximum energies of 
accelerated electrons depend not only on the distance from 
the base of the jet but also on the phase of the binary sys- 



tem. They increase with the distance from the base of the jet 
as oc z 1 / 2 , when the acceleration process is saturated by syn- 
chrotron energy losses (see also Bosch-Ramon et al. 2006). 
However, at the regions of the jet where the magnetic field 
is relatively weak (farther from the base of the jet but still 
close to the massive star), the acceleration process is sat- 
urated by IC energy losses. Then, the maximum Lorentz 
factors of electrons are significantly lower (note the broad 
parabolic dips in Fig. 1, especially evident during the peri- 
astron passage of the compact object). 

In order to achieve the maximum Lorentz factors, ob- 
tained from the above formulae, the acceleration process has 
to occur locally in the jet, i.e. the advection time scale of 
relativistic electrons moving with the jet plasma has to be 
longer than the acceleration time to the maximum energies 
determined by energy losses. The advection (escape) time 
scale of electrons can be estimated by, 

Tadv ~ Z/Vj » 50z/Vj S, (10) 

where the velocity of the jet i)j is in units of the velocity of 
the light, and the distance along the jet z is in units of the 
stellar radius r*. The advection time scale introduces the 
upper limit on the Lorentz factor of the electrons under the 
condition of their local acceleration, 

7max = Accr adv /(mC 2 ) « 1.8 X lO^/Vj. (11) 

For the distances from the base of the jet down to z = 100, 
and acceleration efficiencies in the range £ = 0.03 — 0.3, the 
acceleration of electrons occurs locally since 7m'ix is larger 
than the maximum Lorentz factors given by Eqs. |7| |H| [5] 
(shown in Figs. 0. Regarding particle cooling, it is dom- 
inated by radiative processes for the relevant electron en- 
ergies, providing the particle local cooling, i.e. close to the 
acceleration place. 

By comparing the rate of energy losses of electrons in 
the synchrotron process and the IC process in the Thomson 
regime (Eqs. [3] and |KJ , we find out that at distances from 
the base of the jet lower than zs/t ~ 0.57r]r/T£ , the syn- 
chrotron energy losses dominate over the IC energy losses. 
Therefore, essential contribution to the IC spectra comes 
only from the part of the jet above ~ O.lrv since zs/t is 
typically lower for the range of considered parameters. The 
part of the jet above z m i n is considered since losses prevent 
particles from reaching TeV energies for z < z m in. We fix 
this lower bound on r m i n = O.lr* (for arguments see e.g. 
BOG). From another site, the upper bound on this region, 
Zmax, is introduced by the condition of advection timescales 
being shorter than radiation timescales. Based on Eq. ??, in 
all numerical calculations we fix the upper bound of the jet 
on z ma x = 10r*. Note moreover, that electrons might be in- 
jected by the acceleration mechanism only at specific range 
of distances from the base of the jet. However this possibil- 
ity can not be reasonably specified at the present state of 
knowledge. 

For the acceleration model defined above, we consider 
two scenarios for the injection rate of relativistic electrons 
into the jet. The spectrum of injected electrons is assumed 
to be of the power law type described by the single spectral 
index independent on the distance from the base of the jet z. 
The injection rate (number of particles at specific distance 
z) is assumed to be 
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Figure 2. Spectral Energy Distributions (SED) from cascades initiated by primary electrons injected in the jet with: the efficiency 
depending on the distance from the base of the jet as N e (z) oc 1/z 2 (independently on the phase of the binary system), and differential 
power law spectrum of electrons with the index equal to -2. The upper figures show the 7-ray spectra, produced at the periastron passage 
of the compact object for the observer located at the inclination angles i = 30° and 60°, from the jet moving towards the hemisphere 
containing the observer (the jet) and the sum of the cascade spectra produced by electrons accelerated in the jet and in the counter-jet 
(direction opposite the jet). The jets propagate perpendicular to the plane of the binary system. The 7-ray spectra produced at the 
apastron passage of the compact object are shown in the bottom figures. The specific spectra are calculated for the acceleration conditions 
in the jet described by the parameters: £ = 0.03 and r\ = 0.1 (full curves), £ = 0.3 and r\ = 0.01, (dashed), and £ = 0.3 and rj = 0.1 
(dot-dashed). 



(i) constant along the jet from its base at z m in, 

(ii) drops along the jet as oc z~ 2 . 

We have shown above that the cooling time scale of 
electrons which are able to produce 7-ray photons with en- 
ergies above 1 GeV is shorter than the advection time scale 
of electrons in the jet (see Eq. (13) and the text below). 
Since this cooling time scale of electrons is also much shorter 
than the characteristic time scale of the considered picture 
(determined by the part of the orbital period of the binary 
system), we can assume that all the energy transfered to 
electrons from the acceleration mechanism at specific time 
is radiated almost during the same time. The spectrum of 
7-rays can be calculated correctly by complete cooling of the 
injected electron spectrum. 

The model with continuous injection rate of electrons 
produce very flat 7-ray spectra, inconsistent with the upper 
limits of the Whipple telescope for LSI +61° 303 (Hall et al. 
2003, Fegan 2005)). Therefore, we show only the results of 
calculations for the second injection model of electrons into 
the jet. 

Since the optical depths for very high energy 7-rays in 
the radiation field of the massive star in LSI +61° 303 are 
large, the 7-ray spectra, emerging from the binary system 
toward the observer, are formed in the IC e ± pair cascades 
occurring in the anisotropic radiation of the massive stars. 
We make use of the Monte Carlo cascade code described 
in detail in previous works (Bednarek 2000, B06). However 
the code has been modified in order to include the effects 
connected with the synchrotron cooling of electrons injected 
into the magnetic field of the jet. Therefore, in this work we 
are also able to show the simultaneous synchrotron photon 
spectra produced by electrons accelerated in the jet. When 
calculating the 7-ray spectra escaping from the binary sys- 



tem we neglect the effects of Doppler boosting due to rel- 
ativistic motion of the jets since for the parameters typical 
for LSI +61° 30 (jet speed, ~ 0.6 and observation angles 30° 
and 60°) the Doppler factors are rather small, D m 1.66 and 
1.14, respectively. These relativistic effects change slightly 
the angular and energy distribution of electrons in the jet 
(as observed from the compact object and the observer rest 
frame) and so they are neglected in the present calculations. 
Note moreover that the possible Doppler boosting concerns 
only the primary (first generation) 7-rays produced by elec- 
trons in the jet but not to the secondary cascade 7-rays. 



3 THE CASCADE 7-RAY SPECTRA 

The distance of the compact object from the massive star in 
the binary system LSI +61° 303 changes in the range from 
r p = 1.5r* (periastron) uptor a = 9.15r* (apastron). For the 
periastron passage the optical depths for 7-rays injected at 
most of the hemisphere (on the e ± pair production process 
in 7 + 7 — > e ) are larger than unity (see e.g. calculations 
in Dubus 2006, B06). At the apastron passage the optical 
depths for the 7-rays are larger than unity in some range 
of angles around the massive star (see. e.g. Bednarek B06). 
Note that isotropic electrons in the jet produce 7-rays pref- 
erentially in directions where they find the largest optical 
depth. Therefore, the absorption of 7-rays and the interac- 
tion of secondary e pairs have to be taken into account 
when calculating the 7-ray spectra produced by electrons in 
the radiation field of the massive star both at the perias- 
tron and the apastron passages. We have performed Monte 
Carlo calculations of the 7-ray spectra escaping from the 
binary system at arbitrary directions, applying the above 
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mentioned models for the acceleration of primary electrons 
in the jet, i.e. the injection of electrons along the jet with 
efficiency dropping from the base of the jet according to 
oc z~ 2 . The IC e ± cascade model described in a more de- 
tail in Bednarek (B06) has been applied. In this model we 
assume that secondary cascade leptons are isotropized by 
the random component of the magnetic field in the wind of 
the massive star. Some motivation for that to possibly oc- 
cur is discussed in Bednarek (1997, Sect. 2 of that paper). 
In this cascade process we consider only IC scattering of 
the radiation coming from the massive star neglecting the 
possible contribution from the soft photon field produced in 
the accretion disk. To allow a simple analysis, it is assumed 
that electrons are accelerated with the power law spectrum 
and differential spectral index equal to —2. The injection of 
relativistic electrons in the jet as a function of the distance 
from its base has been also made in other papers concerning 
microquasars (e.g. Bosch- Ramon et al. 2006). The cut-off in 
the electron's spectrum depends on the location of the ac- 
celeration place in the jet and on the phase of the compact 
object. It is determined by the synchrotron, or IC energy 
losses in the jet and by the acceleration efficiency (defined 
by the parameter £). The values of rj and £ are selected 
in order to allow acceleration of electrons up to TeV ener- 
gies (see Fig. 1). Since two jets are expected from a single 
compact object, we consider the production of 7-rays in the 
jet which propagate above the plain of the binary system 
(i.e. the jet directed towards the observer) and in the jet 
propagating below the plane of the binary system (i.e. the 
counter-jet). We note that 7-rays produced in the counter- 
jet may have problems with passing through the matter and 
radiation of the accretion disk around the compact object. 
The transfer of 7-rays through the region of the disk radi- 
ation depends on many unknown accretion disk parameters 
(see e.g. Bednarek 1993) and the column density of mat- 
ter through the accretion disk in LS I 61°+303 is unknown. 
Due to these uncertainties, we show separately 7-ray spec- 
tra from the jet (the case when the counter-jet 7-rays are 
absorbed) and from the jet + counter-jet (the counter-jet 
7-rays are not influenced by the accretion disk). 

In order to obtain the first generation of 7-ray spec- 
tra from the injected electron spectrum we cool injected 
electrons up to energies below 3 GeV applying numerical 
methods. Such approach guaranttee that all the energy of 
electrons is transfered to radiation provided that the obser- 
vation time is much longer than the cooling time of electrons 
below 3 GeV. Due to this complete cooling of electrons in 
the jet, the power in electrons accelerated above 3 GeV has 
be approximately equal to the power in the observed 7-ray 
spectrum above 1 GeV. The results of the calculations of 
the 7-ray spectra from the cascade developed by these first 
generation of 7-rays produced in the jet are shown in Fig. 2 
for the case of an observer located at inclination angles of 
30° and 60°, and the compact object at the periastron (up- 
per figures) and the apastron passages (bottom) (spectrum 
of electrons normalised to 1 MeV and the injection rate to 
one electron between z m i n and z m ax). The 7-ray spectra, 
produced by electrons injected in the jet at the periastron 
passage, steepen significantly at energies above a few tens of 
GeV due to the efficient cascading. This cascading effects are 
smaller when the compact object is close to the apastron. 
The spectra have similar fluxes in the GeV energy range but 
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Figure 3. The multi-wavclcngth spectrum observed from LSI 
+61° 303 (X-ray data from the ROSAT (R - Goldoni & 
Mereghctti 1995) and the OSSE (O - Grove ct al. 1995), and in 7- 
rays from the COMPTEL (C - van Dijk ct al. 1996), the EGRET 
(E - Kniffen et al. 1997), and the detection by the MAGIC tele- 
scope (Albert ct al. 2006)) arc compared on the upper figure with 
the inverse Compton and synchrotron spectra calculated in terms 
of IC e* pair cascade modelling (for £ = 0.1 and r/ = 0.1 and the 
electron spectrum oc E~ 2A ). The IC and synchrotron spectra 
are shown for the periastron passage (dotted curves), and for the 
phase 0.5 at which the maximum 7-ray flux is predicted (dashed 
curves) . The 7-ray light curves (relative units) in the energy range 
1 — 10 GeV and > 200 GeV are shown in the middle figures for the 
inclination angles of the binary system i = 30° (full curve) and 
i = 60° (dashed). The cascade 7-ray spectra escaping towards 
the observer for selected phases of the compact object, these two 
inclination angles and f = 0.1 and r) = 0.1 arc shown in the bot- 
tom figures: the periastron (full curve) and the apastron passages 
(dot-dashed), and the maximum 7-ray flux at phase 0.5 (dashed). 



differ significantly in the TeV energy range for different pa- 
rameters £ and n. Spectra produced at the periastron have 
lower fluxes at TeV energies due to the absorption of 7-rays. 
The contribution of the counter-jet to the total spectrum 
is important only at GeV energies at the periastron by in- 
creasing the flux calculated from the jet by a factor of ~3-4. 
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4 THE CASE OF LS I 61°+303 

In order to be consistent with the observational constraints 
given by the 7-ray flux observed from the EGRET source 
2EG J0241+6119 (Kniffen et al. 1997), and recently mea- 
sured spectrum at TeV energies (Albert et al. 2006), in our 
modelling a steeper spectrum for electrons injected into the 
jet has to be applied (photon differential spectral index equal 
to —2.4). We calculate the 7-ray spectra escaping towards 
the observer for different phases of the compact object in 
order to find out exactly at which phases the highest fluxes 
of GeV and TeV emission should be expected. In Fig. 3 the 
upper panel, the results of cascade calculations (for £ = 0.1 
and 77 = 0.1 and the inclination angle 60°) are compared 
with the multi-wavelength spectrum of LSI +61° 303. The 
IC spectra are shown only above 1 GeV due to the tech- 
nical problems with the calculation of the cascade 7-rays 
spectra at lower energies (restrictions on computing time). 
In fact, the spectra should extend also through the MeV 
energy range (observed by the COMPTEL detector) and 
break at lower energies (in the region of hard X-rays) due to 
the inefficient cooling of electrons in the jet (as discussed in 
Sect. 2). The 7-rays in the GeV- TeV energy range are pro- 
duced in the IC process and the lower energy photons (from 
X-rays up to the soft 7-rays) result from the synchrotron 
emission of electrons accelerated in the jet. We show the 
spectra at the peak of the TeV 7-ray emission occurring 
at phase ~0.5 (dashed curves) and at the periastron pas- 
sage (dotted curves) at which the lowest TeV fluxes are pre- 
dicted, after normalisation to the average GeV flux of the 
EGRET source 2EG J0241+6119 reported in the direction 
of LSI +61° 303. Note that predicted GeV flux could vary 
by only ~ 50%. So such normalization does not have very 
strong consequences when estimating the phase dependent 
TeV 7-ray fluxes. The 7-ray spectrum at TeV energies is 
consistent with the observations by the MAGIC telescope 
at energies above ~200 GeV (the power law spectrum with 
the differential spectral index —2.6 up to a few TeV, Al- 
bert et al. 2006). The MAGIC telescope was not able to 
detect this source at the low emission state (phases close to 
the periastron). According to our calculations this low state 
should be characterised by a flatter spectrum and the flux 
on a level of ~5 lower than at phase ~0.5. Note moreover, 
that the level of synchrotron emission produced at phase 
0.5 is similar to the ROSAT X-ray observations reported by 
Goldoni & Mereghetti (1995) but stronger than the flux re- 
ported from direction of LSI +61° 303 by the OSSE detector 
(Grove et al. 1995). 

In Fig. 3 the middle panel, we show the expected 7-ray 
light curves from LSI +61° 303 in the GeV (1-10 GeV) and 
the TeV energies (> 200 GeV). Clear anti-correlation is pre- 
dicted between the fluxes in these two energy ranges. This 
is due to the importance of the cascading processes inside 
the binary system. The lowest TeV fluxes are expected at 
the periastron passage of the compact object and the highest 
fluxes should be observed after the periastron (with the peak 
at phase ~ 0.5). This specific phase distribution is due to the 
absorption effects of the TeV 7-rays (cascading effects) and 
the location of the observer at the azimuthal angle ~ 70° 
after the periastron. Note that the observer located at the 
inclination angle 30° should observe a high level of TeV emis- 
sion in broader range of phases (half width between ~ 0.4 



up to ~ 0.1), than for the inclination of 30° (between ~ 0.4 
up to ~ 0.9). This feature might serve as a diagnostic of the 
inclination angle of the binary system. The GeV flux reaches 
a peak at phases ~0.15, when the compact object is behind 
the massive star for the observer located at larger inclination 
to the plane of the binary system (i.e. at 60°). For this lo- 
cation, the cascade, developing towards the observer, meets 
the largest optical depths which results in the production 
of large fluxes of the GeV 7-rays. The GeV flux could vary 
by ~50% with the phase of the binary system. Interestingly, 
close to these phase (~0.2) there should also be observed a 
small local peak in the TeV 7-ray light curve. This peak is 
produced by electrons accelerated at larger distances from 
the base of the jet at which the escape conditions of 7-rays 
to the observer are more favorable. Then, TeV 7-rays are di- 
rected at relatively small angles respect to the massive star 
photons which enhances the efficiency of 7-ray production. 

In Fig. 3 the bottom panel, we also show the 7-ray spec- 
tra for a few selected phases of the compact object and both 
inclination angles i = 30° and 60°. It is clear that the GeV 
fluxes and spectra do not vary strongly but the TeV fluxes 
and spectra (> 200 GeV) change significantly with the phase 
of the compact object. Our calculations indicate that the 
TeV spectrum at the periastron should be a factor of ~5 
lower than at the phase ~0.5. Moreover, the TeV spectrum 
at the periastron should be steeper below ~200-300 GeV 
and significantly flatter above this energy respect to the TeV 
spectrum expected at the phase ~0.5. 



5 DISCUSSION AND CONCLUSION 

The massive binary systems provide unique conditions for 
investigating of the basic problem of astrophysics, i.e. the 
acceleration process of particles, due to the presence of a 
well-defined target (the radiation field of the massive star). 
However in order to derive useful constraints the details of 
the sometimes complicated radiation processes by acceler- 
ated particles (e.g. cascades in the anisotropic radiation of 
the massive star) have to be considered. Calculations of the 
type described in this paper predict observational features of 
the 7-ray emission (in the GeV to TeV energy range) which, 
if observed, strongly indicate on a very efficient acceleration 
of electrons already close to the base of the microquasar jets 
and allows to constrain the efficiency of this acceleration 
process. In fact, the 7-ray emission from the microquasar 
type massive binary LS 5039 has been already observed in 
the TeV energies (Aharonian et al. 2005) and in the GeV 
energies (Paredes et al. 2000), and from LSI +61° 303 is 
observed in the TeV energies (Albert et al. 2006) and in the 
GeV energies (Kniffen et al. 1997). Our cascade IC pair 
model predicts the 7-ray light curve in the GeV energy range 
from LSI +61° 303 (Fig. 3) which shows maximum close to 
the periastron as suggested by Massi (2004). In contrast, the 
TeV 7-ray light curve peaks at the phase ~ 0.5 showing clear 
anti-correlation with the GeV emission (see also earlier cal- 
culations of the propagation of 7-rays inside LSI +61° 303 
by Bednarek (B06)). This cascade model also predicts the 
flattening of the 7-ray emission above a few hundred GeV 
and lower 7-ray flux at phases close to the periastron. 

The 7-ray spectra shown in Figs. 2 and 3 are calcu- 
lated only from the part of the jet extending below lOr* (in 
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Figure 4. The comparison of the 7-ray spectrum produced by 
electrons in the jet directly to the observer (dotted curve) with the 
7-ray spectrum escaping to the observer in the cascade process 
considered in this paper (full curve). The parameters of the exam- 
ple model are the following: £ = 0.1 and r\ = 0.1 and the electron 
spectrum <x E~ 2 , inclination of the binary system i = 60° and 
the compact object is at the apastron passage. 

order to fulfil the condition of local, efficient cooling of elec- 
trons) and above O.lrv Note, however, that the considered 
region of the jet is broad enough as to provide the dominant 
contribution to the escaping 7-rays (see also Bosch-Ramon 
et al. 2006). The region below O.lr* might contribute only 
to the GeV part of the spectrum. Moreover, the IC model 
in which electrons in the inner part of the jet scatter soft 
radiation coming from the accretion disk should predict dif- 
ferent 7-ray light curves. For example, in our model the 
7-ray flux depends on the relative orientation of the jet, 
massive star and the observer predicting the peak of TeV 
emission at phase ~0.5. These features can be tested by the 
future GeV and TeV 7-ray observations. As we mentioned 
in Sect. 1, the IC model in which collimated electrons scat- 
ter stellar radiation has been already considered in a few 
papers, e.g. recently by Bosch-Ramon et al. (2006) , Der- 
mer & Bottcher (2006). These previous works do not take 
into account the IC e pair cascade effects which, as we have 
show above, can play an important role. In order to quantify 
this effect we compare in Fig. 4 the 7-ray spectra produced 
by electrons in the jet directly toward the observer and the 
7-ray spectra which are formed in the cascade process for 
the apastron passage where differences between these two 
spectra are expected to be the smallest. Although the 7-ray 
spectra above a few hundred GeV are almost this same, they 
differ significantly at lower energies. Therefore, the broad 
band shape of the 7-ray spectra are clearly different. Note 
that similar effects are expected for the inclination angle 
i — 30° since differences between the spectra at the apastron 
passage for these two inclination angles are small (see the 
light curves for these two inclination angles in Fig. 3 for the 
phase of 0.73). The cascade process considered in this paper 
significantly re-distribute the energy in the 7-ray spectrum 
toward the observer, influencing e.g. the normalization of the 
spectrum at lower energies. In the works mentioned above 
the contribution of 7-ray emission from electrons scattering 
of the disk radiation in the jets of microquasars have also 
been included. It was shown that this contribution can not 
dominate over the 7-rays produced by scattering the stellar 
radiation. 

On the other hand, electrons accelerated in the jet above 
lOr* might contribute to the highest energy part of the 7- 



ray spectrum and to the low energies (X-rays) as a result of 
synchrotron process and synchrotron self-Compton mecha- 
nism since synchrotron energy losses of electrons might dom- 
inate far away from the massive star (e.g. Atoyan & Aharo- 
nian 1999). However, significant variability of the TeV flux 
with timescales similar to the period of the binary system 
has been recently detected by the MAGIC telescope (Albert 
et al. 2006). So then, the contribution to the observed 7-ray 
spectrum from the outer parts of the jet should be rather mi- 
nor in case particle injection were constant along the orbit. 
The problem of the phase dependent accretion onto the com- 
pact object has not been considered in the present work (for 
possible indications how it can occur see e.g. Bosch-Ramon 
et al. 2006, Christiansen et al. 2006). We concentrated only 
on the radiation processes occurring during uniform injec- 
tion of energy along the jet (independent on the phase). 

The electrons cooled up to a few GeV in the inner jet 
are advected to the outer parts of the jet. They can be re- 
sponsible for the lower energy emission observed from some 
microquasars. Based on the considered here cascade model, 
we estimate that a significant fraction of the particle energy 
(~ 10 ) survives the 0.1 — lOr* region to allow the forma- 
tion of radio jets, but it's treatment is out of the scope of 
the paper. 
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